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ABSTRACT 

We analyse those objects in the Brightest 55 sample of clusters of galaxies which have a 
short central cooling time and a central temperature drop. Such clusters are likely to require 
some form of heating. Where clear radio bubbles are observed in these clusters, their energy 
injection is compared to the X-ray cooling rate. Of the 20 clusters requiring heating, at least 
14 have clear bubbles, implying a duty cycle for the bubbling activity of at least 70 per cent. 
The average distance out to which the bubbles can offset the X-ray cooling, rh ca t, is given 
by rhoatAcooi = 0.86 ± 0.11 where r coo i is defined as the radius as which the radiative 
cooling time is 3 Gyr. 10 out of 16 clusters have rheat Acooi ^ 1, but there is a large range 
in values. The clusters which require heating but show no clear bubbles were combined with 
those clusters which have a radio core to form a second sub-sample. Using rh ca t = 0.86r coo i 
we calculate the size of an average bubble expected in these clusters. In five cases (3C129.1, 
A2063, A2204, A31 12 and A3391) the radio morphology is bi-lobed and its extent similar to 
the expected bubble sizes. A comparison between the actual bubble size and the maximum 
expected if they were to offset the X-ray cooling exactly, i2 max , shows a peak at -Rbubble ~ 
0.7-Rmax with a tail extending to larger -Rbubbie/-Rmax- The offset from the expected value 
of i?bubbic ~ Kmax may indicate the presence of a non-thermal component in the innermost 
ICM of most clusters, with a pressure comparable to the thermal pressure. 

Key words: galaxies: clusters: general - galaxies: clusters: cooling flows - X-rays: galaxies: 
clusters 
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1 INTRODUCTION 

Since the discovery of "hol es" in the X-ray emis : 
tre of the Perseus Cluster l B ohringer et alj[l993lv 
the launch of Chandra, many more depressions in the intra-cluster 
medium (ICM) of low redshift clusters have been found (e.g. Hy- 
dra A.lMcNamara et all200cl A20 52.lBlanton et all200lt A 2199. 
llohnston^UuTl27)02T ~Centaurus. |SandenH^abja al2002l. Re- 
cent c ompilation s are given in[5un^^JJ^200^lDum^^Fabian 
120041) : iBirzan et alj |2004). Such holes have been observed to 
anti-correlate spectacularly with the radio emission from the active 
galactic nucleus (AGN) at the centres of these clusters. Their mor- 
phology, particularly in the closest clusters, has led to the interpre- 
tation that these are bubbles of relativistic gas blown by the AGN 
into the thermal ICM. This relativistic gas is less dense than the 
ICM and so the bubbles are expected to detach fr om the core and 
rise up buoyantly through the cluster, e.g. Perseus IChurazov et alJ 
l200ri iFabian et alj|2003bl) . The older, detached, bubbles tend not 
to have GHz radio emission associated with them and have been 
termed "Ghost" bubbles. 

The X-ray emission of the ICM naturally leads to the con- 
clusion that the plasma should be cooling. To maintain pressure 
support, the gas is expected to flow on to the central galaxy 
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as a "cooling flow." The subsequent increase in density would 
lead to a "cooling-catastrophe", with extremely rapid cooling in 
the cluster centre. However, with the high spatial and spectral 
resolution of Chandra and XMM-Newton little of the expected 
X-r ay emitting cool gas has been found IPeterson et alJ 120031 
see IPeterson & Fabian! 120061 for a review). Many mechanisms 
have been proposed by which the cool gas could be heat ed, in- 
cluding, for example , thermal conduction iKim & NaravanlfeoOSt 
IVoigt & Fabiarl2004h . but this appears not to work for clusters be- 
low 5 keV. 

A majo rity (71 per c ent) of "cooling core" clusters harbour ra- 
dio sources (Burns 199C), and a similar fraction of clusters which 
require heating (likely to be a cooling core) harbour clear bubbles 
(Dunn et al. 2005). The action of creating the bubbles at the cen- 
tre of the cluster by the AGN is a favoured method of injecting 
energy into the central regions of the cluster and so prevent the 
ICM from cooling. This process sets up sound/pressure waves in 
the ICM , the dissipation of which requires the ICM to be viscous 
jFabian et all2003tl l2006 ) . The viscous dissipation of the pressure 
waves allows the energy from the bubble creation to be dissipated 
far from the cluste r centre, as gentle, continuous and distributed 
heating as required l Voigt & Fabian 2004). Birzan et al. 1 2004) sur- 
veyed 80 clusters in the Chandra archive without discriminating 
between cooling and non-cooling clusters and found 16 that con- 
tained bubbles. This was interpreted as clusters having AGN/bubble 
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Table 1. Cluster Sample 



Cluster 


Redshift 


Obsld 


Exposure (ks) 


UjJJTV/u 


0349 


919 


19 7 


3C 129 1 


0223 


2219 


9.6 


A85 


0521 


904 


38 4 


A262 


0164 


2215 


28 7 


A478 


0882 


1669 


42 4 


A496 


0223 


93 1 


18 9 


A 1644 


0474 


2206 


18 7 


A 1 65 1 


086 


4185 


9.6 


A 1689 


181 


5004 


19 9 


A 1795 


0627 


493 


19.6 


A2029 


0767 


89 1 


19 8 


A2052 


0348 


890 


36 8 


A2063 


0355 


6263 


16 8 


A2199 


0.0300 


497 


19.5 


A2204 


0.1523 


499 


10.1 


A2597 


0824 


922 


39.4 


A31 12 


0746 


2516 


16 9 


A3391 


0545 


4943 


18 4 


A3558 


0.0475 


1646 


14.4 


A4059 


0478 


897 


40 7 


AWM7 


0.0172 


908 


47.9 


Centaurus 


0.0109 


4954 


89.1 


Cygnus A 


0.0570 


360 


34.7 


Hydra A 


0.0522 


4969 


96.9 


Klem44 


0.0283 


4992 


33.5 


M87 


0.0037 


2707 


98.7 


MKW3s 


0.0449 


900 


57.3 


Ophiuchus 


0.028 


3200 


50.5 


Perseus 


0.0183 


3209 


95.8 


PKS 0745-191 


0.1028 


508 


28.0 



All the clusters in the sample, for the sample selection see text. The expo- 
sure time is that after reprocessing the data. 

activity only 20 per cent of the time. In this work we take a sam- 
ple of clusters which forms a subset of an almost complete sample 
and investigate the bubble population therein. As the parent sample 
is almost complete we are able to investigate prevalence of cool- 
ing core clusters, duty cycles of bubble activity and radio sources 
within the general cluster population. 

The sample selection is described in Section [2] and the data 
preparation and reduction in Section[3] The subset is itself then split 
into two: those clusters which are expected to require some form 
of heating and harbour clear bubbles; and a combination of those 
which require some form of heating but do not harbour bubbles and 
those clusters which have a central radio source. The analysis of 
those clusters with clear bubbles and of those without is discussed 
in Sections [4] and [5] respectively. The overall implications of this 
work is discussed in Section |6| with future investigations outlined 
in Section^ 



imaged by previous X-ray satellites, and is nearly flux complete at 
high galactic latitudes (all but nine are at b > j20'circ|). 

Out of the B55 sample, two clusters did not have ROSAT ob- 
servations, and four only had ROSAT HR1 observations. Out of the 
remainding 49 clusters, iDunn et alj J2005) selected those which 
had a short central cooling time (< 3 Gyr) from the ROSAT PSPC 
and a large central temperature drop (T out or/7ccntrc > 2), as these 
would be the clusters requiring some form of heating so that large 
quantities of cool gas are not formed. Of the 49 clusters, 23 have a 
short t coo \ and, using the most recent temperature profiles, 21 have 
a large central temperature drop. 20 clusters have both, of which 14 
(70 per cent) show clear depressions (bubbles) in the X-ray emis- 
sion from the Intra-Cluster Medium (ICM). The duty cycle of bub- 
bling in clusters which require heating is therefore also at least 70 
per cent. From the whole sample we find that the total fraction of 
clusters with clear bubbles from current observations is at least 25 
percent (14/55). 

Of the six clusters which require heating but have no clear 
bubbl es, only one h as no radio source at the centre - AWM7 
iFurusho et al 120041) . 2A 0335+096 has a complicated core struc- 
ture, but there are depressions which have been interpreted as bub- 
bles iMazzotta et al 120 03) and so this cluster is included with those 
which have clear bubbles. MKW3s also has a faint depression fur- 
ther out from the cen tre, which has been interpreted as a bubble 
( Mazz otta et alj|2002l) and this too is included with the clusters 
which have bubbles. 

Also taken from the B55 sample are those clusters which are 
not already in the sample and which host a radio core at the centre, 
regardless of whether a ROSAT observation exists. The NVSS 1 was 
used to find central radio sources. The NVSS has extensive, but 
not complete sky coverage, and is not very deep, so clusters which 
harbour faint central sources will be missed 2 . Out of the clusters 
not already in the sample, at least another 10 which have a central 
radio source in the NVSS. So the fraction of clusters which have a 
central radio source is at least ~ 53 per cent (29/55, i.e. including 
those with bubbles). These two cluster samples contain a total of 
30 clusters (see Table 0, 16 of which have a radio source, clear 
bubbles and require heating 3 , 13 which have a radio source but no 
bubbles, and AWM7 which has no evidence for radio activity. The 
assignment of the clusters to the different groups is shown in Table 

HI 

We have chosen a central cooling time of t coo i < 3 Gyr to in- 
vestigate those clusters where a cooling flow would form if there is 
no heating. These are therefore the clusters in which extreme heat- 
ing rates are expected. Changing to t coo i < 7 Gyr ~ tHubbie/2 
could add another nine clusters into the sample which require heat- 
ing. There is only one cluster of these nine which also has a central 
temperature drop - A2142 - and so the fraction of clusters with 
bubbles from those which require heating changes little Three of 
these nine clusters (A2063, A1689 & A1651) already fall into the 
sample which have core radio detections. None has a report of clear 
bubbles. 



2 SAMPLE SELECTION 

We use an upd ated version of the sample first outlined in 
IDunn et at]l2005t) . which is a subset o f the Brightest 55 (B55) sam- 
ple. The B55 sample was compiled by Edge et al. 1 1990) and stud- 
ied in detail with ROSAT data from a cooling-time point of view by 
IPeres et alj Jl998l) . The B55 sample is a 2 - 10 keV flux-limited 
sample of clusters of galaxies which are close enough to have been 



1 NRAO (National Radio Astronomy Observatory) VLA (Very Large Ar- 
ray) Sky Survey 

2 For example. Ilvlarkovic et al. 12004) present new detections of central 
sources in A1650 and A2142 which are not detected in the NVSS. 

3 iDonahue etal] l2005t) investigate both A 1650 and A2244 as they have 
radio-quiet coolin g cores (t r . o i < H^ 1 ). A2244 has t coo \ > 3 Gyr from 
IPeres et at] (l99§) and athough Mark ovic et alj l2004tl detect a radio source 
at the centre of A1650 none is seen in the NVSS and so neither is in our 
sample. 



© 0000 RAS, MNRAS 000, 000-000 



Investigating AGN Heating in a Sample of Nearby Clusters 3 



Table 2. Cluster Sub-samples 



Clear Bubbles NVSS Radio No Radio 



Heating Heating No Heating Heating 



A85 2A0335+096 1 " 3C 129.1 AWM7 

A262 A496 A 1644 

A478 A2204 A1651 

A1795 PKS 0745-191 A1689 

A2029 MKW3s+ A2063 

A2052 A3 112 

A2199 A3391 

A2597 A3558 

A4059 Klem44 

Centaurus Ophiuchus 

Cygnus A 

Hydra A 

M87 

Perseus 

14 5 10 1 



NOTES These clusters have complicated cores/faint depressions and are 
added to the Clear Bubbles set (for further discussion see text). The Clear 
Bubbles (incl t ) clusters form one subset and the remainder form the other. 
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2 4 6 8 10 12 14 16 18 



U,i (Gyr) 

Figure 1. The distribution of t coo i in the B55 sample. The black bars in- 
dicate those clusters with clear bubbles, the grey are for the clusters with 
no bubbles but with radio sources. This graph has been truncated, two clus- 
ters with t coo i = 22 and 32 Gyr (A1367 & A2255) have been removed to 
expand the rr-axis. 

2.1 Heating Distribution 

The distribution of the central cooling times of the B55 clusters is 
shown in Fig. Q Almost all of the clusters which have a cooling 
time less than 3 Gyr contain a radio source or clear bubbles, with 
the exception being AWM7. The clusters in the B55 sample which 
do not have central cooling time values from the ROSAT PSPC have 
not been included in the figure. 



3 DATA PREPARATION 

The X-ray Chandra data of the clusters were processed and cleaned 
using the CI AO software and calibration files (CIAO v3.3, CALDB 
v3.2). We began the reprocessing by removing the afterglow detec- 



tion and re-identifying the hot pixels and cosmic ray afterglows, 
followed by the tool ACIS_PROCESS_EVENTS to remove the pixel 
randomisation and to flag potential background events for data ob- 
served in Very Faint (VF) mode. The Charge-Transfer Efficiency 
was corrected for, followed by standard grade selection. Point- 
sources were identified using the WAVDETECT wavelet-transform 
procedure. For clusters observed with the ACIS-S3 chip, the ACIS- 
Sl chip was used to form the light curves where possible. In all 
other cases, light-curves were taken from on-chip regions as free 
as possible from cluster emission. For the spectral analysis, back- 
grounds were taken from the CALDB blank-field data-sets. They 
had the same reprocessing applied, and were reprojected to the cor- 
rect orientation. 

Cluster centroids were chosen to lie on the peak in the X-ray 
surface brightness. Annular regions were automatically assigned 
with constant signal-to-noise, stopping where the calculated sur- 
face brightness of the cluster dropped below zero. The initial signal- 
to-noise was 100, and this was increased or decreased by successive 
factors of \/2 to obtain a number of regions between four and ten. 
The minimum signal-to-noise allowed was 10. 

The 0.5 — 7keV spectra were extracted, binned with a 
mini mum of 20 c ounts per bin, and, using XSPEC (vl2.2.1r) 
(e.g. |A rnaud 1996J), a PROJCT single temperature MEKAL (e.g. 
iMewe et al] 1 19951) model with a PHABS absorption was used to 
deproject the cluster. In some clusters the temperatures for some 
of the regions were undefined. To solve this problem, the mini- 
mum number of regions was reduced, a maximum radius for the 
outermost annulus was set, or the outermost region was removed 
to try to improve the behaviour of the profile (this was required 
for 3C 129.1, A496, A 1689, A2063, A2597, A3391, A4059, Klem 
44, Ophiuchus & PKS0745). Using the deprojected cluster tem- 
perature, abundance and normalisation profiles; density, pressure, 
entropy, cooling time and heating profiles were created. The tem- 
perature profiles and the derived profiles for the clusters with clear 
bubbles are shown in the Appendix (Figs.l9lto ll5t . 

These profiles give aziumthally averaged values for the clus- 
ter properties and have been used in the subsequent calculations. 
In some clusters, notably M87 and Perseus, the central parts of 
the clu ster are not very smooth, e.g. due to bubbles. [Donahue et alj 
feOOfjt) show that these features do not strongly bias estimates of the 
entropy. As such the use of these azimuthally averaged values is not 
likely to introduce large biases into the subsequent calculations. 



4 CLUSTERS WITH BUBBLES 

To estimate the energy input of the AGN to the ICM we estimated 
the energy required to create the observed bubbles. If the expansion 
rates are slow then this is the sum of the bubble's internal energy 
and the PthdV work done. 

E= — — -Pt h V + P th V= — ^—PthV (1) 
71 — I 71 — I 

where V is the volume of the bubble, P t h is the thermal pressure of 
the surrounding ICM and 71 is the mean adiabatic index of the fluid 
in the bubble. In the case where the fluid in the bubble is relativistic, 
71 = 4/3, and so the total energy in the bubble is 4P t hV. As the 
bubbles are approximately elliptical we parameterise with a semi- 
major axis along the jet direction, Hi, and a semi-major axis across 
it i? w . The volumes are therefore 4ttR\R^/3, where we have as- 
sumed that the bubbles are prolate ellipsoids. 

To obtain an estimate for the rate at which energy from a bub- 
ble would be dissipated in the ICM, the age for the bubble, t ag c, 
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is required. Most of the bubbles in the sample are young and still 
attached to their radio core. In these cases the timescale used is the 
sound speed timescale (i agc = t Cs = R/c s ) where the sound speed 
is given by 

c s = \J -y 2 kBT/^niH (2) 

where 72 = 5/3 for a non-relativistic gas. This timescale arises 
from the fact that no strong shocks have been observed in clusters. 
From this the bubbles are assumed to be expanding at not much 
more than the local sound speed of the ICM. In the cases where 
the bubbles are ghost bubbles, i.e. they do not contain GHz radio 
emission and are detached from the cluster centre, we have used the 
buoyancy timescale (t agc = tb = R/vb) as these are expected to 
be rising buoyantly up through the ICM, with the buoyancy velocity 
given by 

v b = y/2gV/SC D (3) 

where Cd = 0.75 is the drag coefficient IChurazov et alJl200lh 
and S = nrl, is the cross-sectional area of the bubble. Those bub- 
bles which have been classed as Ghost bubbles are indicated in 
Table [3] For furt her di s cussion on the t i mesca les see Section l6~2l 
and lDunn et alJ (2005); Dunn & Fabian 12004). The energies and 
ti mescales for the bubbles were calculate d from values for the radi i 
in lAUenetall <2006l) ; lDunn et alJ <2005t): jDunn & Fabianl {2004 . 
except for 2A 0335+096 which are from birzan et alJ feOOfl) ! 

From the density profiles and the relevant cooling function, 
the heat input required per spherical shell was calculated. By com- 
paring this to the total power provided by the bubble (Pbubbic = 
4PthVyt a g C ), the radius out to which the bubble power can offset 
the X-ray cooling was calculated for each cluster (rh ca t, see Ta- 
ble|5J- The radius is quoted in kpc and also as a fraction of r coo i 
(for a cooling time of 3 Gyr). For an estimate on the uncertainties, 
Monte-Carlo simulations of the calculations were performed. 

The bubbles in A2052, Cygnus A and Perseus supply enough 
energy to offset all the X-ray cooling far beyond the cooling radius. 
In fact the energy they supply offsets all the X-ray cooling in the 
analysed regions of the cluster. In the Perseus cluster this is mainly 
because the cluster is so close, and so only the innermost regions 
fit onto one chip. In A2052 and Cygnus A the surface brightness 
of the X-ray emission falls to zero before the edge of the chip. As 
a result the spectral analysis stops at a radius such that the bubble 
energy offsets more than all the X-ray cooling within that radius. 
In the case of M87, the cluster is so close that the cooling radius, 
r coo i, could not be calculated from a single chip; however the bub- 
ble energy is greater than the X-ray cooling within the largest radius 
obtained on a single chip. These clusters are highlighte d in T able 
The results for M87 obtained from iGhizzardi et alJ <2004 are 
shown for comparison and indicate that the current jet and counter 
jet cavities approximately offset all the X-ray cooling within r coo i. 
This result, however, has not been included in any of the subsequent 
analysis except where explicitly stated. 

The average distance, as a fraction of the cooling radius, out 
to which the bubbles can offset the X-ray cooling is estimated us- 
ing the mean of the results in Table [3] The uncertainties in the 
mean were estimated using a simple bootstrapping method. As a 
fraction of the cooling radius the mean distance out to which the 
bubbles offset the X-ray cooling is ri lcat /r coo i = 0.86 ± 0.11. 
The mean power that the bubbles supply, as a fraction of the X-ray 
cooling within the cooling radius, is 0.89 ± 0.16. The distribution 
of the distance out to which the bubbles offset the X-ray cooling is 
shown in Fig. [3] The clusters where the bubble power is such that 
it is greater than the all X-ray cooling within the analysed region 




Lcooi (erg/s) 



Figure 2. The X-ray Luminosity within the cooling radius versus the bub- 
ble power for the clusters. Clusters that fall above the line have AGN which 
inject sufficient energy into the ICM to offset all of the X-ray cooling. The 
four clusters with open circles are those with Ghost bubbles (2A 0335, 
A25 97, A85, A478 and MK W 3s). The gTey point is M87 using the data 
from lGhizzardi et alJ 120041) . 



(A2052, Cygnus A & Perseus) and M87 have not been included in 
the estimate on the mean. Adding A2052, Cygnus A and Perseus 
into the estimation of the mean gives rv at /r coo i = 1.45 ± 0.29 
(with a heated fraction of 3.49 ± 1.93), and so it is likely that the 
value of ri lca t /r coo i = 0.86 is a lower limit for this sample of clus- 
ters. See Section|6|for further discussion of the range of rh ca t /r coo \ 
and its implications. Comparing L C ooi with Pbubbie (Fig.|2j shows 
that there is a slight correlation between the cooling luminosity and 
the bubble power. The clusters lying above have AGN which supply 

sufficient energy into the ICM to offset the cooling. 

In a set of clusters harbo uring clear bubbles iBirzan et alJ 

120041) and lRaffertvetalJ 12006) found that the mechanical lumi- 
nosities required to offset the cooling ranged between lP t hV and 
20Pth V. Around half the objects in their sample had cavities which 
(assuming 4P t h V) could offset the cooling though there is a large 
spread, with some objects falling short. In this sample, which is 
drawn from an almost flux-complete parent sample, we find similar 
results. 

As yet the method by which the bubbles dissipate their en- 
ergy into the ICM has not been clearly determined. The energy also 
needs to be transported out to a large radius, for example, the vis- 
cous dissipation of sound waves in the ICM as they travel out in 

■ III 1 _| 

the cluster 1 Fabian et al. 2003a, 2005). As these are seen at large 

radii, they can plausibly carry their energy sufficiently far out. As 
a result, the bubbles analysed here, on average, supply enough en- 
ergy to offset the X-ray cooling, but whether this energy actually 
performs this function is currently unknown. 
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Table 3. Bubble Heating 



Cluster 


Bubble 


Ghost? 


n 




^dist 


L coo? 


r ool 


fheat 


pat. / 

" ' cool 


Heated Fraction 2 








(kpc) 


(kpc) 


(kpc) 


10 43 ergs -1 


(kpc) 


(kpc) 




2A 0335 


E 


G 


9.3 


6.3 


23 


22.7 ±0.4 


69.8 ±1.0 


14.8 ± 1.62 


0.21 ±0.02 


0.09 ±0.02 




W 


G 


4.8 


2.6 


28 












A85 


N 


G 


5.3 


7.0 


14 


9.0 ±0.4 


47.6 ±1.3 


81.4 ± 10.7 


1.71 ± 0.23 


2.54 ±0.71 




S 


G 


6.4 


8.6 


22 












A262 


E 




2.4 


1.8 


3.2 


0.5 ±0.0 


31.4 ±0.8 


28.9 ± 3.20 


0.92 ± 0.10 


0.93 ±0.10 




W 




2.7 


1.9 


3.4 












A478 


NE 


G 


4.0 


6.9 


9.8 


78.2 ± 3.0 


84.8 ±2.7 


18.8 ± 2.63 


0.22 ± 0.03 


0.13 ±0.02 




SW 


G 


4.5 


7.7 


10.7 












A1795 


NW 




4.1 


3.1 


3.8 


26.9 ± 2.1 


68.8 ±3.3 


38.1 ± 1.63 


0.55 ± 0.04 


0.35 ±0.04 




S 




4.8 


2.8 


5.2 












A2029 


NW 




7.2 


2.2 


9.4 


51.2 ± 9.4 


73.3 ±9.9 


31.9 ± 3.50 


0.43 ± 0.08 


0.31 ±0.08 




SE 




6.5 


2.9 


9.4 












A2052 


N 




10 


10 


10 


4.5 ±0.1 


46.2 ±1.0 


226 ± 0.00 


4.90 ±0.11 


6.83 ±0.83 




S 




15 


15 


15 












A2199 


E 




2.3 


1.3 


2.3 


5.3 ±0.3 


46.0 ±1.1 


16.7 ±0.92 


0.36 ± 0.02 


0.20 ±0.03 




W 




2.4 


1.4 


2.4 












A2597 


NE 


G 


7.8 


7.8 


21 


33.0 ± 1.5 


66.6 ±1.0 


63.2 ± 5.08 


0.95 ± 0.08 


0.85 ±0.14 




SW 


G 


12 


7.8 


25 












A4059 


N 




5.4 


2.2 


5.4 


2.4 ±0.2 


45.5 ±2.0 


46.91 ±3.22 


1.03 ± 0.08 


1.06 ±0.18 




S 




3.9 


4.8 


3.9 












Centaurus 


E 




2.8 


1.7 


2.5 


2.1 ± 0.0 


42.8 ±0.2 


20.7 ±2.99 


0.48 ± 0.07 


0.44 ±0.07 




W 




2 7 


2 


3 












Cygnus A 


E 




30 


21 


44 


29.8 ± 1.1 


55.1 ± 1.6 


359 ± 0.00 


6.52 ± 0.19 


38.0 ±5.8 




W 




34 


22 


48 












Hydra d 


N 




20 


5.5 


18 


15.6 ± 0.2 


56.2 ±0.7 


83.5 ± 11.6 


1.49 ± 0.21 


1.42 ±0.14 




S 




17 


4.4 


22 












M87 


E 




1.5 


1.5 


1.5 


2.5 ±0.0 


27.4 ±0.0 


27.4 ±0.00 


1.00 ±0.00 


0.66 ±0.02 




W 




1.9 


1.2 


1.9 












MKW3s 


S 


G 


13 


21 


56 


5.0 ±0.3 


53.4 ±2.0 


110 ±29.7 


2.06 ± 0.56 


2.48 ±0.67 


Perseus 


N 




8.2 


8.2 


8.2 


50.2 ± 0.1 


76.4 ±0.1 


109 ± 0.00 


1.43 ± 0.00 


2.29 ±0.03 




S 




8.9 


8.9 


8.9 












M87 e 


E 




1.5 


1.5 


1.5 


1.5 ±0.1 


37.4 ±2.3 


35.3 ±3.0 


0.94 ±0.10 


0.93 ±0.15 




W 




1.9 


1.2 


1.9 













The clusters in bold are those where the bubbles heat beyond the analysed region. For M87 the cooling radius could not be determined from a ACIS single 
chip. These clusters were not included in the calculation of the mean values. 

a L coo \ calculated for the 0.5 — 7.0 keV range. b The radius, as a fraction of the cooling radius, out to which the bubble power (4P t j 1 V/t a ge) can offset the 
X-ray radio source by interacting with its surroundings cooling. c The fraction of the X-ray cooling that occurs within the cooling radius which is offset by the 
energy of the bubble's expansion. d The bubble sizes for Hydra A a re not those from the cl uster scale outburst from Nulsen et al. 1 2005) but from the smaller, 
inner radio emission. e These values for M87 come from the data in iGhizzardi etafl <2004l) which extend to larger radii. 



5 CLUSTERS WITHOUT BUBBLES 

The second cluster subset is a combination of those clusters which 
require heating and show no clear bubbles, and those clusters which 
harbour a radio source (see Table [2j. The amount of X-ray cool- 
ing within ri lcat /r coo i = 0.86 (the mean distance out to which 
the bubbles offset the X-ray cooling for the other cluster subset) 
for these clusters was calculated. In the assumption that any X-ray 
cooling from these clusters is also offset by AGN bubbles, and that 
these bubbles are young and so expanding at the sound speed, the 
dimensions of these bubbles was calculated. The resulting dimen- 
sions are shown in Table [4] where we have assumed that bubbles 
occur in pairs. The uncertainties have been estimated using simple 
Monte-Carlo simulations of the calculation. 

The total X-ray counts expected within the area of the pre- 
dicted bubble in the current observation was calculated from the 
average X-ray counts per pixel in the central regions of the clus- 
ter. The central regions were chosen so that it could be reasonable 
assumed that the surface brightness is constant. This is, however, 



a lower limit on the number of counts and the signal to noise for 
the innermost regions, i.e. those where any bubbles are likely to be. 
This gave the expected X-ray signal to noise which is shown in Ta- 
ble|4] Using the bubbles in A2052 and Hydra A as templates, we 
estimate that the counts in the centre of the bubbles are around 30 
percent lower than the counts in the rims. 

So that any bubbles could be detected above the noise at a 3<r 
level the X-ray signal to noise has to be greater than around 10 
(equivalent to noise at 10 per cent level). This does, however, rely 
on knowing where the bubbles are in the centre of the cluster. If 
there is extended radio emission then the average counts interior 
and exterior to the expected bubble can be compared to see if there 
is a net decrement within the bubble. Half of the clusters in the 
sample have signal-to-noise values of less than 10, and so, even 
if extended radio emission is detected, no significant detection of 
bubbles from the current X-ray observations are likely. The counts 
per pixel in the central regions of these clusters is also very low. 

A496, A2204, A3 112, AWM7, Klem44, Ophiuchus and 
PKS 0745 all have X-ray signal-to-noise which is such that if a 
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rheat/l"cool 

Figure 3. The distribution of the radius (as a fraction of the cooling radius) 
out to which the energy in the bubbles can offset the X-ray cooling. The 
extra points at r hcat /r cool = 0.8-1 and r hcat /r cool = 1.4 - 1.6 
are M87 and Perseus respectively. A2052 and Cygnus A are not included 
on this plot. All four of these clusters were not taken into account when 
calculating the mean, indicated by the vertical line. 



bubble of the expected size were present it should be detectable in 
the current X-ray observation. This is, however, still with the caveat 
that the size and location of the bubble is known (again, extended 
radio emission would help with this). 

To be able to detect a decrement in the X-ray emission without 
any radio emission to guide the eye is more difficult. AWM7 and 
Klem44 only have around one count per pixel in the central regions, 
and so it is unlikely that a bubble could be identified by eye from 
the X-rays alone. It is more likely that in clusters with higher counts 
per pixel (A2204, Ophiuchus and PKS 0745) that bubbles would be 
detectable as these have more than 1000 counts expected within the 
area of the expected bubble. 

The contrast ratio we have chosen, however, is for very clear 
bubbles with well defined rims in nearby clusters. In smaller ra- 
dio sources, bubbles may be less clearly defined, and so the X-ray 
signal-to-noise may have to be much higher to be able to clearly 
detect cavities if they are there. 

Using the Very Large Array (VLA) 4 and Australia Telescope 
Compact Array (ATCA) 5 archives we attempted to find observa- 
tions of these clusters to investigate the morphology of the radio 
sources at their centres. In most cases the radio sources were not 
resolved in the archival observations. However for at least four of 
the clusters in Table|4]there are extended radio sources which could 
be bubbles that have not been detected in the X-ray images of clus- 
ters. In other cases it has not been possible to find radio data of high 



4 The National Radio Astronomy Observatory is operated by Associated 
Universities, Inc., under cooperative agreement with the National Science 
Foundation. 

5 The Compact Array is part of the Australia Telescope, which is funded by 
the Commonwealth of Australia for operation as a National Facility man- 
aged by CSIRO. 
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Figure 4. The 5 GHz radio emission from 3C 129.1 [Taylor et all200lli . 



enough resolution, and so there may be further clusters with, as yet, 
undetected bubbles. 

In the case of 3C 129.1 the morphology of the radio emis- 
sion is very suggestive that the surrounding ICM is constraining 
the expansion of the inner radio lobes (see Fig.|4j- The X-ray ob- 
servations to date, however, show no clear interactions of the radio 
source with the ICM. The observed radio lobes in 3C 129.1 are 
~ 2.6 kpc (5.5arcsec) in radius, which is around four times that 
expected for an "average" bubble in this cluster, 0.60 ± 0.16 kpc 
(1.34 ± 0.36arcsec). The observed radio lobes correspond to 
around 10 Chandra pixels in radius (which are 0.49 x 0.49 arcsec), 
and so if the radio lobes are excavating cavities in the ICM, then the 
resultant bubbles should be seen in deep X-ray observations of the 
cluster which would increase the X-ray signal-to-noise. 

A3391 also harbours an extended radio source w ith a bi- 
lobed morphology lOtani et afll9 98 ; Morga nti et al J 19991) . The ra- 
dio lobes in this cluster are separate from the radio core. If these 
lobes have formed cavities in the ICM, this separation implies 
that the cavities are older bubbles which have detached from the 
cluster centre and are rising buoyantly. The radio structure is not 
smooth and simp le, but the 13 cm (2.38 GHz) observations from 
lOtani et alj \ 19981) show the lobes as ~ 8 x 12 kpc (8 x 12 arcsec 
or 16 x 24 Chandra pixels) in size. The expected size of an average 
bubble is 2.61±0.33 kpc, which is much smaller than the observed 
size of the radio emission. 

The ROSAT X-ray maps of the centre of this cluster show no 
clear indication that the radio source is interacting with its sur- 
roundings. An explanation for these large detached radio lobes is if 
close feedback exists between the X-ray cooling and the AGN heat- 
ing, then there may recently have been a past cycle of AGN heating 
and now there is a period of quiescence and no (small) lobes are 
seen with the expected sizes. 

The central radio source in A2063 also has two peaks of 
emission. The dimensions of the observed radio emission are ~ 
2x1 arcsec corresponding to ~ 1.4 x 0.7 kpc (or 4 x 2 Chandra 
pixels). The expected sizes for average radio bubbles in this cluster 
are 2.09 ± 0.38 kpc (2.95 ± 0.53 arcsec or ~ 6 Chandra pixels) 
in size. 
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Table 4. Predicted Bubble Sizes 
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^*cool 


Predicted Bubble Radius 


X-ray S/N 


counts/pix 




(kpc) 


(kpc) 


arcsec 
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r: qq _|_ o 


0.60 ±0.16 


1.34 ±0.36 
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A496 


54 QQ _|_ 1 Q7 


5.50 ±0.53 


8.40 ± 0.80 


53.7 


3.15 


A 1644 




1.98 ±0.34 


2.13 ±0.36 


4.6 


0.37 


A1651 


1 4 24 + 1 2 53 


1 AC\ -I- AQ 


u.o / m u.ou 


1.5 


0.22 


A 1689 


63.12 ±5.85 


5.19 ±0.42 


1.70 ±0.14 


7.0 


1.32 


A2063 


29.27 ± 1.57 


2.09 ±0.38 


2.95 ±0.53 


6.5 


0.37 


A2204 


69.44 ±3.16 


8.13 ± 1.04 


3.06 ±0.39 


34.1 


9.56 


A3112 


57.86 ±2.71 


5.44 ±0.62 


4.05 ± 0.46 


30.4 


4.34 


A3391 


10.84 ±2.18 


2.61 ±0.33 


2.52 ±0.32 


2.5 


0.08 


A3558 


14.90 ±1.46 


1.60 ±0.10 


1.72 ±0.10 


4.0 


0.41 


AWM7 


24.84 ± 2.34 


1.25 ±0.22 


3.58 ±0.63 


15.6 


1.46 


Klem44 


38.79 ± 1.41 


2.56 ±0.21 


4.61 ±0.39 


16.1 


0.94 


Ophi 


32.85 ±0.82 


1.69 ±0.32 


3.01 ±0.57 


28.6 
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90.76 ±2.16 
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Figure 5. The 4.9 GHz radio emission from A2063. 



Figure 6. The 1.5 GHz radio emission from A2204 ISanders et all2005t) . 



None of the X-ray observations of these three clusters are deep 
enough so that clear indications of an interaction between the radio 
source and the ICM are expected. 

The 1.32 GHz radio map of A3 112 shows two faint diffuse 
region s of emission t o the south east and south west of the core 
( Takiz awa et aljl2003i) . Unfortunately the radio image has an elo- 
gated beam which makes it difficult to determine whether a faint 
excess in the X-rays is definitely associated with the radio lobes. 
The observed radio emission is ~ 4 arcsec in radius correspond- 
ing to ~ 5 kpc (or ~ 8 Chandra pixels). The emission is offset 
from the core, so it may correspond to detached bubbles which still 
have some GHz radio emission. The expected sizes for average ra- 
dio bubbles in this cluster match almost exactly the observed radio 
emission. 

A2204 and PKS 0745 are the two of the more distant clus- 
ters in our sample, with redshifts greater than 0.1 (0.1523 and 
0.1028 respectively). The radio morphology of A2204 shows two 
clear maxima (see Fig.|6j, which could be the two lobes of the ra- 



dio source, which have not been fully resolved. Again, the X-ray 
emission shows no clear indication of any interaction of the radio 
source with the ICM. The radii of the observed radio maxima in 
A2204 are ~ 5 kpc (~ 2 arcsec or 4 Chandra pixels), and the size 
of the average bubble expected in this cluster is only slightly larger 
(8.13 ± 1.04 kpc or 3.06 ± 0.39 arcsec). 

PKS 0745 has an amorphous radio source at the cluster centre, 
rather than a clear bi-lobed morphology and it is unclear as to the 
effect of the central radio source on the ICM. The total dimensions 
of the observed radio source are 9.2 x 6.0 kpc, which is not that 
dissimilar from that for an expected bubble (r — 10.1 ± 2.0 kpc). 

The X-ray observations pf these three clusters are sufficiently 
deep that if any bubbles existed and the surface brightness contrast 
between the rims and the bubble centre were only 10 per cent they 
should be observed, especially as there is extended radio emission 
to guide the eye. It is possible that the contrast between the bubbles 
and the surrounding rims in these two clusters is very low, and as a 
result any cavities have not yet been detected. In any case, deeper 
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X-ray observations would allow further investigation into the pres- 
ence or absence of bubbles in these clusters. 

The current X-ray observations for most of these clusters are 
too short to provide sufficient signal-to-noise for any cavities to 
be clearly identified in the images. As a result some cluster may 
harbour bubbles which cannot be detected. Some of the clusters do, 
however, harbour extended radio sources whose dimensions are not 
dissimilar to those expected from an "average" bubble, should one 
be present in the cluster. 



6 DISCUSSION 

We first discuss the spread in values of r / r coo i and then investigate 
the possibility of the young bubbles being below their maximum 
size. Finally we detail some of the uncertainties and biases present 
in this analysis. 

Although the mean of the radius out to which the bubbles can 
offset the X-ray cooling is 0.86r coo i (from bootstrapping), there is 
a large spread in the values (see Fig.|5J. In some clusters the bubble 
power offsets the heating out to over 1.5r coo i (Hydra) whereas in 
others it only reaches out to 0.36r coo i (A2199). The range of the 
energy supplied by the bubbles as a fraction of that required within 
r coo i also has a large spread, from xl.4 that required in Hydra, 
to xO.20 in A2199, with an average of 0.89. The Ghost bubbles 
are not taken into account here, along with the clusters where the 
analysis is not complete (A2052, Cygnus A and Perseus, but not 
M87). Adding these in the range increases from 0.22 to 6.5r coo i. 
So in some clusters the bubbles can provide (more than) enough 
heating to offset the X-ray cooling, but some fall short. 

Fig. Q shows the distance to which the bubbles can offset the 
X-ray cooling versus the central cooling time (from IPeres et all 
11993 adjusted to our cosmology). T he gre y point is M87 using 
the data presented in Ghizzardi et al. (2004). The three squares are 
A2052, Cygnus A and Perseus where the bubbles provide energy 
to beyond where it was possible to deproject the cluster. Unfortu- 
nately there is no clear correlation between r/r coo i and t CO oi- 

Seven clusters have r/r coo i > 1, three have r/r coo \ ~ 1 and 
six have r/r coo i < 1. So in most (10 out of 16) of the clusters 
the energy supplied by the bubbles is sufficient to offset the X- 
ray cooling. Of the six which have insufficient, four have young, 
active bubbles. A possible explanation is that these bubbles are still 
growing and so have not reached their maximum size. Therefore 
they currently do not contain sufficient energy to offset the cooling, 
but may do at later stages in their evolution. For further discussion 
see Section l6~T1 

The ghost bubbles are more difficult to explain as these bub- 
bles would be expected to be at their final size. They are also those 
which supply the least amount of energy into the ICM. The cool- 
ing times at the centres of these clusters are very short (< 1 Gyr) 
and so a comparatively large amount of heating is required to offset 
the X-ray cooling. It is possible that we have caught the AGN in 
a period of relative quiesence where the ICM is cooling and only 
small bubbles are produced. Once a reservoir of fuel for the AGN 
has built up from the ICM then if it goes into a major outburst the 
bubbles may be such that they offset the cooling. 



6.1 Bubble Sizes 

So far we have assumed that the bubbles in the clusters are at their 
"maximum" size, however is this true? Calculating the bubble ra- 
dius expected if its energy were offsetting all the X-ray cooling we 
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Figure 7. The distance out to which the bubbles can offset the X-ray cool- 
ing versus the central cooling time. The three clusters which are have lower 
limits are Perseus, A2052 and Cygnus A (in increasing t 0DD i). The four 
clusters with open circles are those with Ghost bubbles (2A 0335, A2597, 
A85, A478 and MKW 3s). The gTey point is M87. 



show that it is likely that our sample contains a large fraction of 
small bubbles. 

If the jet supplying the bubbles is constant, then there is ex- 
pected to be a continuous cycle of bubble growth with periodic de- 
tachment. As the radius of the bubble depends on the cube-root 
of the bubble age, most bubbles wo uld be observed when t hey are 
close to their maximum sizes. From Chura zov et alJ 12000). in the 
assumption that the jet is of constant power and the bubbles are 
expanding subsonically, 
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CLt 

7i 
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3CL 



167rP t hV 



rav 



(4) 



We naively assume that the numerical constant, C = 1, 71 = 4/3 
and t max = -Rmax/tJgrav Using the sound speed, v Cs , for w grav as 
this depends only on the surrounding ICM and not the bubble size, 
and half the power required to balance cooling for L — L coo \/2 
(as there are two bubbles per cluster), P ma x is calculated. This is 
the maximum size of the bubble possible if all the energy of its 
creation is converted to PthV work on the ICM. The distribution of 
Pav/Pmax is shown in Fig. [8] where P av = \/R\R'i,, the average 
radius of the observed bubble. Monte Carlo simulations of the data 
have been used to take into account the errors in the distribution 
and provide an estimate on the values of the bins 6 . All clusters 



6 The uncertainties on ij a v / Rmax result from the uncertainties in the in- 
put data, e.g. temperature profiles of the cluster. To obtain an estimate on 
the uncertainties on the values of the binned J? a v / Pmax distribution, the 
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Figure 8. TOP: The distribution of bubbles as a fraction of their maxi- 
mum size (i?/fl max ) assuming subsonic expansion. BOTTOM: The result 
of Monte Carlo simulations of the distribution in order to obtain errors (the 
vertical lines) on the distribution, and so a clearer indication of the true dis- 
tribution. The lighter grey bars are from M87. The bubbles from Cygnus A 
and one from A2052 would appear off to the right of this plot. 



with young bubbles (11 clusters, 22 bubbles) are included in the 
analysis and Fig. [8] As our analysis of M87 did not go far enough 
out in radius to ob tain r coo i we used the temperature, and electron 
density profiles in Ghizzardi et al. 1 2004) to calculate i? av / i? max . 
The values of t max and -R av /-R ma x are shown in Table|5| 

The distribution peaks around ~ 0.7i? av /i? max with a tail 



results were binned many times (10 4 runs), assuming a Gaussian distribu- 
tion of uncertainties in ij av / -R ma x ■ The inter-quartile range in the values 
for each of the bins is what is shown by the error bars. 



extending to larger i? av /i? ma x- This implies that a large number of 
the bubbles have not yet attained their maximum size, and others 
are much larger than they "should" be. 

The bubbles with a small _R av /-Rmax are likely to be young 
bubbles caught in the act of growing. They are not particularly 
younger than any other bubbles, but may be young for the clus- 
ter. For example, the inner bubbles in A2199, which have been 
analysed here, are very much smaller than the large scale radio 
emission, which also corresponds to clear decrements in the X-ray 
emission from the ICM. The large X-ray holes have i? av ~ 15 kpc 
giving an fl av /-Rmai ~ 5, whereas the inner bubbles have R av ~ 
1.6 kpc giving an -R av /-Rmax ~ 0.5. These bubbles are also those 
which do not supply sufficient energy to the ICM, which is unsur- 
prising as we have used L coo i when calculating ii max . 

As the bubbles with i? av /-Rmax < 1 are likely to be young 
bubbles, it is possible that those bubbles with i? av /-R ma x > 1 may 
be those which have detached and are expanding as they rise up 
through the ICM. 

Also, some of the clusters with i? av /-R ma x > 1 have mor- 
phologies which are different to the younger ones. Hydra A may not 
have "bubbles" in the sa me way as the other clusters. The cluster- 
scale outburst reported in Nulsen et al ] <2005l) shows that the radio 
source has ploughed through the ICM rather than being confined 
by it; Cygnus A and A4059 are similar - there are X-ray decre- 
ments, but the synchrotron plasma is not totally confined by the 
ICM. It may be the case that the bubbles in A2052 and Perseus are 
just in the process of detaching, they still appear as young bubbles 
as the radio emission has not had time to decay, but they are at the 
maximum possible size for the cluster and are beginning to buoy- 
antly rise. The bubble dimensions we have used for the Perseus 
cl uster in this analys is are not those labelled as "Inner Bubbles" 
in lDunn et alj l2W 6). Using the radii for the Inner bubbles in this 
analysis (~ 4 kpc) gives an -R av /-Rmax ~ 0.6. 

If we assume that most bubbles that we observe are at their 
maximum size, then we have overestimated i? max . To move the 
peak in the distribution, such that i? av /-R ma x ~ 1, Bmax has to 
reduce by a factor ~ 0.7. Which of the assumptions that we have 
used in the analysis could be changed to reduce 7? max ? 

We have assumed that the timescale over which the bubbles 
heat the ICM is the same as the sound speed timescal e of the bub- 
bles' formation. Bubbles may grow in fits and starts iFabian et alJ 
2005), the actual time taken to expand the bubble may be shorter, 
as there are periods where there is no expansion. Depending on the 
exact nature of the growth the true expansion time may be less than 
the one calculated. The energy dissipation time is also going to be 
different to the bubble creation timescale (see Section l6"2l for more 
discussion on this topic). 

We have assumed 71 =4/3, but if it is possible for 71 ~ 1.1 
then the peak of the distribution occurs at -R av /-Rmax ~ 1. A 
smaller value of 71 means that the same amount of energy creates a 
smaller bubble than if 71 = 4/3. Investigations into the weak shock 
around the northern Perseus bubbles have shown it to be isother- 
mal iFabian et alj|200f3) . which requires an effective 72 ~ 1 for 
the ICM. Further investigations are required to determine the true 
natu re of 7 in the radio bu bbles and the ICM. 

ISanders et all <2005h find non-thermal X-rays from the cen- 
tral regions of the Perseus Cluster. These most likely arise from 
magnetic fields and cosmic rays in the ICM. Using the more recent 
900 ks observation Sanders et al. (in prep) find that the non-thermal 
electron pressure within a radius of 40 kpc is comparable to the to- 
tal thermal pressure. As a result the total pressure against which 
the bubbles are expanding is about twice that which has previously 
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Table 5. r a v/-Rmax 
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8.89 ±0.00 


306 ±2 


8.98 ±0.05 


0. 


65 


± 





.01 


1. 


49 


±0. 


01 


M87 
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The average bubble radius (R av ), ICM pressure, sound speed (V s ) and -Rav/^max for all the young bubbles in the sample. The starred clusters are those in 
which the bubbles do not provide sufficient energy to offset the X-ray cooling. 



been estimated. Amending Equation[4]results in 

^max = \ Tc Ton \ = ?K~ = 0.7i?max- (5) 

y 167r(2P t h)^grav V2 

If the non-thermal components are taken into account, then f? ma x 
reduces by the factor of 0.7 required. As the distribution of 
Rav/Rmax peaks at around 0.7 then it is likely that a non-thermal 
component of the ICM exists and is important in The inner- 
most refions of all cool core clusters. We note that the possible 
presen ce of a n on-thermal pressure component was suggested by 
IVoigt&Fabianl<2006h from determination of mass profiles in cool 
core clusters. 



6.2 Timescales 

The timescales used to calculate the power output into the cluster 
by the creation of bubbles have been those of the sound speed. So 
far no strong shocks have been observed around the rims of the 
bubbles. This means that they are currently expanding at less than 
the sound speed. We therefore have an upper limit on the bubble 
expansion time, which has been used as the timescale over which 
the bubbles deposit energy into the ICM. 

This is the bubble expansion timescale, however, rather than 
the timescale over which the bubbles dissipate their energy in 
the ICM. We have so far assumed that these two timescales are 
equivalent. But, as the bubbles can still be seen in the form of 
detached ghost bubbles (e.g. Perseus & Centaurus), bubbles last 
longer than this timescale. As such the energy is dissipated over a 
longer timescale than has been assumed here. 

If bubbles are continually produced, so the moment one de- 
taches, or shortly thereafter, another one starts to form, then the 



difference in timescales is not a problem. Although the energy dis- 
sipation occurs over a longer timescale than the bubble creation, 
the net input of energy into the central regions of the cluster oc- 
curs on the bubble creation timescale. However, if the duty cycle of 
the central engine is short, so that bubbles are produced only rarely, 
then the energy dissipation timescale is the more relevant. Accurate 
estimates of the time-averaged energy dissipation rate are vital to 
the study of the energy input into the central regions of these clus- 
ters. Our study above indicated that the duty cycle is at least 70 per 
cent in clusters which require some form of heating. As there are a 
number of clusters which, although they do not have clear bubbles, 
have complicated central morphologies or radio emission, the duty 
cycle may be as high as 90 per cent. 

In some clusters, Hydra A for example iNulsen et af]|2005l) . 
there appears to have been a "cluster scale outburst." In these cases 
the outburst may have heated the central parts of the cluster suffi- 
ciently so that it is not, currently, cooling very fast. Therefore the 
observed bubble size and the current cooling rate may not be as 
closely coupled as we have assumed. The heating effect of such 
large scale requires further investigation, as in these cases the ICM 
does not confine the bubble and the ICM pressure cannot be taken 
as uniform. 

There exist radio sources where the morphology indicates 
that a ctivity has restarted in the centre after a time of quies- 
cence l Schoenm akers et alj|200fj : ISaripalli e t al. 2002). Although 
for some of the clusters presented here even a duty cycle of 100 per- 
cent would be insufficient to provide the energy required to off-set 
the X-ray cooling, for others a duty cycle of less than 100 per-cent 
would suffice. 

For further improvements to the line of investigation we have 
taken, the rate of the energy dissipation into the ICM, rather than 
the energy injection, needs to be determined. 
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6.3 Bubble Visibility 

The clusters with short cooling times are also likely to be the ones 
which are X-ray bright as the gas, and hence surface brightness, 
densities are high in the centre. Therefore the bubble-ICM contrast 
is also likely to be high. This could cause a selection effect for 
the percentage of clusters with short central cooling times in which 
clear bubbles have been detected. There may be bubbles in clusters 
without short central cooling times, but as the X-ray emission in 
the centre is not as bright, they have not been detected. 



7 FURTHER WORK 

During this analysis we have assumed that the energy contained 
within the bubble and which is available to the ICM is APthV, as- 
suming that the bubbles are filled with a relativistic gas. Apart from 
energy los t to sound waves, whi ch has been neglected in these cal- 
culations IChurazov et"ail 2002). 4-PthV^ is the maximum amount 
of energy available to offset the X-ray cooling occurring in the clus- 
ter. In some cases this energy suffices to offset the cooling out to 
the cooling radius, in others it falls short. 

If, however, the energy were only P t hV, then in many more 
clusters the AGN would appear not to be able to offset the X-ray 
cooling. As a result, an important question that requires answering 
is whether there is AP^V available in the bubble and whether it 
can all be transferred to/dissipated in the ICM? 

As discussed in Section lo"2l the timescales used here to ob- 
tain the bubble powers are those for the bubble age and not those 
over which the bubbles would liberate their energy into the ICM. 
They also do not measure the time over which sound waves would 
dissipate. These timescales also need to be accurately measured to 
subsequently be able to accurately determine the heating rate of the 
AGN. 

Recent work by Fabia rTet alj 120061) investigating the weak 
shock surrounding the bubbles in the Perseus cluster showed that 
the energy content of the post shock gas was around 2P t hV . The 
observation that the gas across the shock is isothermal also raises 
the possibility that the bubbles could expand faster than the sound 
speed. Thermal conduction across the shock would even out any 
temperature differences, so making the shock less detectable. A 
final issue is the energy in the compressed cosmic rays and mag- 
netic field present in the centres of clusters. Sanders et al. in prep 
show that the thermal and non-thermal pressures are comparable in 
the central regions of the Perseus cluster. As a result, the pressure 
against which the bubbles are expanding may be around double 
what has previously been estimated. The E a ,/i? ma x calculation 
(Section f6.1> implies that this non-thermal component is important 
in other clusters. As these inferences come from a single cluster 
(albeit the X-ray brightest and best studied), further investigation 
is required to determine whether these properties occur in all clus- 
ters and the effect that they have on the assumptions used in AGN 
heating arguments. 



8 CONCLUSIONS 

We have investigated a sub-sample of the Brightest 55 cluster sam- 
ple. At least 36 per cent of clusters have a cool core, and of these 
at least 70 per cent harbour clear bubbles. First we analysed those 
clusters which have a short (< 3 Gyr) central cooling time and a 
central temperature drop of a factor of two. These clusters are ex- 
pected to have the most rapid cooling in the centre and so require 



some form of heating. At least 36 per cent of clusters have a cool 
core, and of these at least 70 per cent harbour clear bubbles, imply- 
ing a duty cycle also of at least 70 per cent. Up to 90 percent could 
have some form of bubbles. The energy provided in the bubbles 
was compared to the X-ray cooling within the cooling radius (for 
r coo i at t coo \ = 3 Gyr). The mean distance out to which the bubble 
power can offset the X-ray cooling is rhcat / r C ooi 0.86 ±0.11, 
with a large spread. In most clusters (10 out of 16) the AGNs en- 
ergy input to the central regions of the cluster is sufficient to offset 
the X-ray cooling. The bubble energy in the remainder does not 
offset most of the X-ray cooling. In some cases this is the result of 
catching these bubbles at a very young age (e.g. A2199). 

Also analysed were those clusters in the B55 sample which 
have a central radio source. The fraction of clusters which have a 
central radio source is at least 50 per cent. These were combined 
with the clusters from the previous subset which did not have clear 
bubbles. For clusters without clear bubbles we have used the mean 
value of rhcat / r C ooi to estimate bubble sizes from the X-ray cool- 
ing. In five clusters (3C 129.1, A2063, A2204, A3312 & A3391) the 
radio sources have bi-lobed morphologies, whose sizes are similar 
from the expected sizes. Only in two of these five clusters is the 
X-ray observation of sufficient depth that any interaction between 
the radio and ICM could be expected to be seen. It is not clear why 
no clear bubbles are seen, especially in A2204. As a result bubbles 
may be more common than previously thought. 

We investigated the ratio of the actual bubble size to that of 
their maximum size if they were to offset all the X-ray cooling, 
and find the peak occurs at 7i av /ii max ~ 0.7 with a tail extending 
to higher Rsv/Rwm- Either we are biased to imaging young bub- 
bles or some assumption is incorrect. A likely explanation is that 
a non-thermal component comparable to the thermal component is 
present in the innermost regions of most clusters. 
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APPENDIX 

Here we present the temperature, density, pressure, entropy, cooling 
time and cumulative heating profiles for the clusters in our sample 
which harbour clear bubbles. 

As can be seen from the figures, the temperature profiles (Fig. 
|9j do not lie on top of each other, but the central temperature drops 
can be seen in these clusters. The electron density profiles are much 
more uniform and much smoother. These features carry into the 
derived profiles; the pressure's do not lie on top of each other and 
fall as expected with radius. 

The entropy 7 and cooling time profiles for the clusters are very 
similar. The scatter about the general trends are much less in these 
profiles than in th e others. The ove r all shape of these profiles are 
similar to those in lPratt et ail feOOfj) : ! Voigt & Fabianl <2004 . 

We have fitted simple poweiiaws to all of the profiles of the 
clusters individually. The mean, maximum and minimum are given 
in the figure captions. In some cases the values for the outermost 
shell have been igno red when thi s is obviously different from a 
global trend I Johnstone et al. 2005). If a shell has not been used for 
a fit in the denisty or temperature, it has not been used any of the 
fits in subsequent derived profiles. 



7 —2/3 

We use S = kTn c as a proxy for the true entropy. 
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Figure 9. The temperature profiles of all the clusters in the sample which 
harbour clear bubbles. Error bars have been omitted for clarity. The 1 — a 
errorbars are ~ 2 — 10 per cent at lOkpc and ~ 0.6 — 15 per cent at 
100 kpc. The average powerlaw index is 0.32 with a range of 0.11 — 0.59. 



Figure 11. The pressure profiles of all the clusters in the sample which 
harbour clear bubbles. The legend is missing deliberately, but is the same 
as the other profiles. The average powerlaw index is —0.60 with a range of 
-0.09 to -1.13. 





r(kpc) 



r(kpc) 



Figure 10. The electron density profiles of all the clusters in the sample 
which harbour clear bubbles. Error bars have been omitted for clarity. The 
1 — cr errorbars are ~ 1 — 3 per cent at 10 kpc and ~ 2 per cent at 100 kpc. 
The average powerlaw index is —0.90 with a range of —0.13 to —1.34. 



Figure 12. The entropy profiles of all the clusters in the sample which har- 
bour clear bubbles. The average powerlaw index is 0.94 with a range of 
0.68 - 1.41. 
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Figure 13. The cooling time profiles of all the clusters in the sample which Figure 15. The profiles of all the clusters in the sample which harbour clear 

harbour clear bubbles. The average powerlaw index is 1.19 with a range of bubbles showing the amount of power required to offset the X-ray cooling 

0.83 - 1.76. perkpc. 




Figure 14. The profiles of all the clusters in the sample which harbour clear 
bubbles showing the cumulative amount of power required to offset the X- 
ray cooling. Note that the differential of these curves, the heating radio per 
kpc, is approximately flat (Fig. 1151 . The average powerlaw index is 1.40 
with a range of 0.80 - 2.09. 
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